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Abstract 

, We study the photoproduction of large-pr charged hadrons in ep collisions, both 

j_i ■ for the inclusive case and for the case where a jet in the final state is also measured. 

Our results are obtained by a NLO generator of partonic events. We discuss the 
sensitivity of the cross section to the renormalisation and factorisation scales, and 
to various fragmentation function parametrisations. The possibility to constrain the 
parton densities in the proton and in the photon is assessed. Comparisons are made 
with HI data for inclusive charged hadron production. 



1 Introduction 



The photoproduction of large-pT jets, photons or hadrons are privileged reactions to 
study QCD and to measure the parton distributions in the proton and in the photon. In 
the past, the interest was mainly focused on the production of jets El, 0, and more recently 
on the production of photons M. Particularly interesting are the dijet cross sections |4], |J, 
or the photon-jet cross sections ||, because the measurement of two jets or particles in 
the final state allows one to constrain the incoming parton kinematics and to explore the 
parton distributions in an accurate way. 

In this paper we present results concerning the large-j?T photoproduction of a charged 
hadron and a jet, ep — > /i ± + jet. This reaction offers several interesting features in 
comparison with dijet or photon-jet reactions. With respect to the latter, the hadron-jet 
cross section is much higher ; the observation of a hadron is also easier than that of a 
photon. Compared to the dijet reaction, the hadron-jet one is also easier to measure, this 
fact being particularly true at small px (~ 5-10 GeV) where it is difficult to model the 
underlying event contribution and to unambiguously measure the transverse momentum 
of two jets. It should also be possible to explore a larger rapidity domain for the hadron 
since there is no cone "hitting the edge" of the detector. 

These features are important when one focuses on the measurement of the parton dis- 
tributions in the photon. As is well known, there are two contributions to photoproduction 
cross sections : the so-called direct and resolved contributions. The resolved contribution 
is important at small values of pt and at large positive rapidities, a kinematic domain 
that the hadron-jet reaction should allow one to explore. 

The theoretical description of the hadron-jet cross sections requires the knowledge of 
the fragmentation functions of quarks and gluons into hadrons. This might appear as a 
disadvantage with respect to the jet-jet reaction. However one has to note that hadroni- 
sation corrections are also needed in jet reactions to describe the evolution of partons into 
hadrons and they are not totally under control. Moreover fragmentation functions are 
now well measured in LEP experiments and several recent NLO parametrisations of quark 
and gluon fragmentation functions are available |7|, [8|, ||]. These fragmentation functions 
have been tested in inclusive charged hadron production in pp, jp and 77 collisions by the 
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authors of ref. @ ; they found a good agreement between theory and data, which confirms 
that the fragmentation functions are under control. 

In this paper we concentrate on the hadron-jet physics and assess the possibility 
to constrain the parton distribution functions of the photon and of the proton. In 
jet-jet reactions, it is usual to constrain the momentum of the incoming partons by 
means of the variables x%Cl = (pj^^e^^ 1 + p^ t2 e^ 3 ^ 2 )/2E p '' y and of the variables 
x*l2 = Pt (e^ 6 ' 1 + e^ 6 ' 2 )/2£' p ' 7 , where E p, ' y are the energies of the incoming pro- 
ton respectively photon. (We follow the HERA convention with the proton momentum 
oriented toward the positive z-axis and the photon momentum toward the negative z- 
axis). For the Born contributions, with only two jets in the final state, the variables x p ' 7 
exactly correspond to the longitudinal momentum fractions of the partons in the proton 
and in the photon. When higher order QCD corrections are considered, the variables x p '" / 
do not fix the initial momenta any more, but they still put useful constraints on these 
momenta. The situation is different in the hadron-jet case because the hadron momentum 
Pjr is only a fraction of the momentum of the outgoing parton, such that x p ' 7 (with p 3 ^ ll 
replaced by ptj,) do not even for the Born cross-section correspond to the incoming parton 
momenta anymore. Therefore it is interesting to study the usefulness of such variables 
in the hadron-jet case, and how they can help to measure the proton and photon parton 
distributions. 

The results presented in this paper are based on a NLO Monte Carlo program which 
generates events containing 2 or 3 partons in the final state. One of them fragments into 
a large-pT hadron, and the others are recombined into jets. This NLO event generator 
provides a flexible approach to implement experimental cuts and to calculate a large variety 
of observables. 

The theoretical framework and the Monte Carlo program are presented in section 2. 
In section 3 we compare theoretical predictions for the inclusive (i.e. no jet observed) 
cross section with HI data ||10|| . The emphasis is put on the study of the sensitivity of the 
cross section to the factorisation and renormalisation scales, and to different fragmentation 
function parametrisations. In section 4 we study the hadron-jet cross section and explore 
the distributions da/dx p '^, in particular their sensitivity to the gluon densities in the 
proton and the photon. Section 5 is the conclusion. 
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2 Theoretical framework 



The NLO Monte Carlo program used in this paper has already been described in refs. 
1 11, ^] in which the photoproduction of isolated prompt photons is studied. Therefore 



we only give a few indications on the general structure of the program here and discuss 
the new features specific to the photoproduction of hadrons. 

In photoproduction events, the electron acts like a source of quasi-real photons whose 
spectrum can be described by the Weizsacker- Williams formula 



fjy) = ^±i 1 + ( 1 -y) 2 ln Qla.ii- v) _ 2(1 - y) ■ (1) 

7 2vr [ y mjy 2 y J 

The quasi-real photon then either takes part directly in the hard scattering process, or 

it acts as a composite object, being a source of partons which take part in the hard 

subprocess. The latter mechanism is referred to as resolved process and is parametrised 

by the photon structure functions F a / 1 (x J ,Q 2 ). Thus the distribution of partons in the 

electron is a convolution 

F a/e (x e , M ) = J o d y dx i fjiv) F a/ji x j, M)8{x 1 y - x e ) (2) 
where in the "direct" case F a / 7 (x 7 ,M) = b~ ai 8{l — x 7 ). 

The production of the final hadron h with momentum Ph is described by a fragmen- 
tation function D%(z, Mp) where z is the fraction of the longitudinal momentum of the 
parton a carried away by the hadron h. The production cross section of a large-py hadron 
and a jet is written as a convolution of the distributions of initial partons, the fragmenta- 
tion of the final parton and the hard scattering cross sections 



da ep ^ hjet (P p: P e ,P h ,P jet ) = J2 J dx eJ dx P J dz F a/e (x e ,M) F b/p (x p ,M) 

a,b,c 

xda ab - c > et (x p P p , x e P e , P h jz, P jet ,v, M, M F )D h c (z, M F ) . (3) 



The hard cross sections do ah ^ c ^ et are calculated at the NLO accuracy. They are expansions 
in powers of ol s {\i) 

da jb->cjet = ^ ^cjet + a 2 (/x) d^cjet^ ^ ^ + ^3) (4) 
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dS*™** = a 2 M dd a £oRN et + ^(M) da a H b o Cjet ^, M, M F ) + 0{pft . (5) 

In expressions ([|), (|j) and @, we have explicitly written the dependence on the large 
scales fi, M and Mp. (For simplicity we choose the same factorisation scale M for the 
incoming photon and proton.) The cross section (||) would not depend on these large 
scales if it were calculated to all orders in a s (fj,). But after truncation of the series (Q) and 
(||), the cross section will depend on /u, M and Mp. The //-, M- and Mp - dependence 
of the HO terms partially compensates the scale-dependence of the Born cross sections, 
and in the next sections we shall study the sensitivity of da ep ~* h i et to the renormalisation 
scale (j, and the factorisation scales M and Mp. 

Let us now discuss the various components of formula (|3|). The fragmentation func- 
tions D%(z, Mp) have been accurately measured in LEP experiments and several NLO 
parametrisations of the latter are now available |7| ||, ^] • There exist non-negligible differ- 
ences between the parametrisations of individual D%(z, Mp), especially at large z > 0.8, 
but for cross sections involving sums over the flavours a and over the hadrons h the differ- 
ences are tiny. In photoproduction, the weights of the different flavours which contribute 
to the cross section are identical to those of the e + e~-annihilation reaction (at least for 
the direct contribution). Therefore the differences between the individual contributions 
should be smoothed when summed over to form the cross section. In the next section 
we shall check the sensitivity of the photoproduction cross section to the fragmentation 
functions. 

We now turn to another important component of formula (|3|), the quark and gluon 
distributions in the photon. In this paper we use a new NLO parametrisation of these dis- 
tributions, the AFG02 parametrisation [13]. The AFG02 parametrisation is an evolution 
of the AFG []14|] parametrisation. The new distributions are more flexible : for instance 
the shape and normalisation of the non-perturbative gluon distribution can be modified. 
This fact allows one to study the sensitivity of cross sections to the gluon distribution. 
The normalisation of the non-perturbative quark distributions can also be adjusted. Con- 
trarily to the AFG distributions, the AFG02 distributions also contain the bottom quark 
distribution. In this paper we shall use the default AFG02 parametrisation which is almost 
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identical to the AFG parametrisation ; the only differences are the number of flavours and 
the value of A^^ = 300 MeV (A^r = 200 MeV was used in the AFG parametrisation in 
agreement with the value of determined ten years ago). As a result, due to a faster 
QCD evolution, the AFG02 distributions are slightly higher at small x and lower at large 
x than the AFG distributions. 

For the parton distributions in the proton, we use MRST99 (g j) |22j as default. For the 



hadron-jet cross section, we will also use the new MRST01 |23| ] and CTEQ6M [24] sets for 
comparison. 

We use a strong coupling constant calculated by solving exactly (i.e. without expansion 
in logQ 2 /A 2 ) the two-loop renormalisation group equation. We work with Nf = 4 active 
flavours as default. Using Nf = 5 the total cross section increases by about 5%. 

Expression ([|) is calculated via a MC code which generates 2 — ► 2 and 2 — > 3 parton 
configurations according to weights given by the subprocess cross sections and the dis- 
tribution functions. A phase space slicing method is used to isolate and to analytically 
calculate the soft and collinear singular contributions of the 2^3 cross sections. The 
soft divergences are cancelled by the corresponding divergences contained in the virtual 
corrections to the 2^2 processes (UV divergences are removed by renormalisation in the 



MS scheme). The remaining collinear singularities are factored out and absorbed in the 



distribution and fragmentation functions using the MS scheme. 

This Monte Carlo code, which uses the event generator Bases/Spring [|l5|], is a par- 
tonic event generator : it contains negative weights coming for instance from the virtual 
corrections to the Born cross sections. It is very flexible and allows one to study various 
cross sections involving a large-pr hadron and jets. Experimental cuts are easily taken 
into account, as well as different jet algorithms. In this paper, we use the fc^-algorithm 



| lq| to define the jets. 



3 Large-px hadron inclusive cross section 

In this section we study the photoproduction of large-pr inclusive hadrons. Similar 
studies and comparisons with data have already been performed in several publications M, 



17, |lq|. Therefore we shall not carry out an exhaustive study of this reaction, but we shall 
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concentrate on the sensitivity of the cross section to the factorisation and renormalisation 
scales, and to different parametrisations of the fragmentation functions. 

In Fig. ||, we display a study of the scale dependence. The theoretical curves are 
compared to HI data [1(1 corresponding to the following kinematic conditions. The pa- 
rameters of the Weizsacker- Williams formula are Q^ax = 0-01 GeV 2 and 0.3 < y < 0.7. 
The ep center of mass energy is y/S ep = 300 GeV. The cross section for large-pr charged 
hadrons is measured in the pseudo-rapidity domain \rj\ < 1. The theoretical curves are 
obtained with the BFGW fragmentation functions ||] and we use AgL = 300 MeV. 
We can see from Fig. p] that the data are described fairly well, in particular that no in- 
trinsic hp is needed to describe the data at low pp, contrarily to what has been found in 
the E706 experiment on fixed target inclusive 7r° production [|l^, B(| . 

In Fig. H we show on a linear scale the ratios of the cross sections calculated with 
M = Mp = [i = C pt (C = 0.5,2) to the cross section calculated at C = 1, where 
Pt is the transverse momentum of the final state hadron. Results obtained with the 
KKP H fragmentation function parametrisations are also displayed^. We clearly see a 
strong dependence of the cross section on the scales for pp ^ 7 GeV. Below 7 GeV the 
perturbative calculation is not reliable, the HO corrections being very large. Moreover, 
for the choice C = 1/2 we explore a range of the factorisation scales (M = Mp = pp/2 ~ 
1.5 GeV) which is very far from the kinematic region where the fragmentation functions 
have been constrained. 



* The dip at pr ~ 5 GeV in the KKP ratio may stem from the fact that the charm threshold is at 
Mf = 2m c in these paramtrisations. Therefore there is no charm contribution below pr ~ 6 GeV in the 
case C = 0.5, whereas for C = 1 the charm contributes already at 3 GeV. In the BFGW parametrisations 
the charm threshold is at Mf = m c , such that the charm contributes at pr > 3 GeV even for C = 0.5. 
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Figure 1: da/drjdp^ for 3 different scale choices compared to HI data. 
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Figure 2: da/dpp calculated with scales \i = M = Mp = C pr (C=0.5, 2) normalised to 
da/dpT at n = M = Mp = pp. 
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We also find a strong scale dependence of the rapidity distribution da/dn which is 
measured in the range 3 GeV < pp < 12 GeV. In Table || we give the results of a study 
of Act = Jo dr]jfc- in which we separately vary u and M, keeping fixed Mp = pp. The 
upper number is the total cross section (direct + resolved) and the lower number is the 
ratio r = HO / {Born + HO). The variation of the cross section with C M = \xjpp is very 
strong for small values of (C M ~ 0.5) corresponding to large values of a s (u). There is 
no region where the cross section o"(Cm, C^) is almost independent of Cm and C^, and the 
ratio r is always large. We again conclude that the theoretical predictions are not reliable 
for pp ~ 3 GeV, a p-p-region which gives an important contribution to the cross section 
integrated over 3 GeV < pp < 12 GeV. 





Cm 


0.5 


1 


2 


C M =0.5 


a 


3.941 


4.509 


4.807 


r 


(0.54) 


(048) 


(0.31) 


C M =1 


a 


2.580 


3.079 


3.430 




r 


(0.57) 


(048) 


(041) 


C M =2 


a 


1.855 


2.264 


2.564 


r 


(0.58) 


(0.52) 


(046) 



Table 1: u, M-dependence of the total cross section integrated over < r] < 0.5 and 
3 GeV< pt <12 GeV. The fragmentation scale has been kept fixed to Mp = pp. The bold 
numbers denote the total cross section in nb, the numbers in italic represent the ratio 
r = HO /(Born + HO). One can see that there is no region of stability for a value of pp m 
as low as 3 GeV. 

The situation improves when we study the sensitivity to scale changes for larger values 
of pp. For instance at pp = 7 GeV - although the sensitivity is still large when we vary 
all scales by the same factor C (see Fig. |2|) - we find a region where the cross section as 
a function of M and n has a flat behaviour. We can see from Table ||, where Mp is kept 
fixed to pp/2, that there is a stability region for O.S^C^^O.b and 0.5<Cm<2. Fixing 
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H = OApT and M = 1.5 pr at the saddle point of the surface given in Table we also 
studied the dependence on Mf and found that the cross section varies by ±2% when Mf is 
varied between 0.3 pt and pr- Therefore the scale sensitivity appears to be under control 
for large values of px {pr 7 GeV) . We shall pursue this study in the next section for the 
hadr on-jet cross section. 





Cm 


0.3 


0.5 


1 


2 


0^=0.3 


a 


0.254 


0.243 


0.229 


0.224 


r 


(0.25) 


(0.06) 


(-0.19) 


(-0.40) 




a 


0.220 


0.222 


0.223 


0.228 


r 


(0.38) 


(0.27) 


(0.13) 


(0.03) 


C M =1 


a 


0.176 


0.184 


0.192 


0.200 


r 


(0.45) 


(0.38) 


(0.29) 


(0.23) 



Table 2: fx, M-dependence of the total cross section integrated over \rj\ < 1 at pr = 7 GeV. 
The fragmentation scale has been fixed to Mp = Pt/2- The bold numbers denote the 
total cross section in nb, the numbers in italic represent the ratio r = HO /{Born + HO). 
In this case there is a rather flat region for small values of \x and large values of M, and 
the overall variation of the cross section due to scale changes is much smaller. 

The sensitivity to the fragmentation function parametrisations ||, ||] is less pro- 
nounced than the one due to scale variations. This result is illustrated in Fig. ||. For px > 
7 GeV, there is at most a 20% difference between the various parametrisations, the BFGW 
and KKP parametrisations being quite close to each other. The BFGW set2 parametrisa- 
tion has a slightly higher gluon than the default BFGW for 0.2 < z < 0.6 and a lower gluon 
for z>0.6, due to a higher N g in the input parameter set. Since the HERA kinematics 
do not probe high z values, the difference to the BFGW default set is negligible in our 
case, as can be seen from Fig. ||| However, the sensitivity to the fagmentation functions 
is non-negligible for pp ~ 3 GeV, a range explored by HI. In Fig. |4] we display a compari- 
son between HI data and theoretical results obtained with different sets of fragmentation 
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Figure 3: da/dpx with different fragmentation functions normalised to da/dpx with 
BFGW (setl) fragmentation functions, at the scales fj, = M = Mp = pt and the hadron 
rapidity integrated over the range \rj\ < 1. 



functions. The dispersion of the results is quite large. In this regard it is instructive to 
look at the distribution da/dz, which is not a physical quantity (the momentum of the 
parton " parent" of the hadron h cannot be measured) , but which gives interesting indica- 
tions on the average value {z) in the cross section, and on the variance (z 2 ) — (z) 2 . The 
^-distributions are displayed in Fig. || for the direct and resolved components. One can 
see that two different ranges in z of the fragmentation functions are probed by the direct 
respectively resolved component. The mean value {z) is larger for the resolved contribu- 
tion^. Since the various fragmentation function parametrisations differ mostly at high z, 
a change of parametrisation will have a different effect on the direct and on the resolved 

^This is due to the fact that the cross section for the production of partons behaves like da / dp^ rton ~ 
A/(p^ rton ) n (n > 0), such that, using p^ rton — p^/z, the behaviour of the cross section is da/dz ~ z"~ 2 
convoluted with the fragmentation functions which decrease with increasing z. The value of n is larger for 
the resolved part, leading to a larger value of (z). 
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Figure 4: da/drj with different sets of fragmentation functions compared to HI data. 

contribution. This effect can be seen in Fig. ^| where we observe that the differences be- 
tween the various parametrisations are more pronounced for positive rapidities, a region 
corresponding to a large resolved contribution. We again verify in Fig. |] that the value 
of p™ in =3GeV is too small to allow a reliable prediction, the dispersion of the results for 
the various parametrisations being too large at low pp. 

We further investigated the impact of using an expansion in logQ 2 /A 2 for a s instead 
of the numerical solution of the two-loop renormalisation group equation, and found the 
result being about 10% lower when using the expansion in logQ 2 /A 2 . We also compared 
our result to the one of Kniehl, Kramer and Potter ||. Using the KKP fragmentation 
functions and the expansion in logQ 2 /A 2 for a s , we obtain the result given in || within 
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Figure 5: Distribution of the hadron momentum fraction z for the direct and resolved 
contributions. (3 GeV< p T < 12 GeV, \rj\ < 1, fi = M = M F = p T .) 



the numerical errors. 
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4 Hadron-jet cross section 



In this section we explore the features of the hadron-jet cross section. The input param- 
eters are the same as for the inclusive cross section, and we consider the cross section 
da/drf 1 drf et (where rf 1 is the pseudo-rapidity of the observed charged hadron), integrated 
in the range 7 GeV < p T < 15 GeV and E 3 ^ 1 > 5 GeV. When there are two jets in the final 
state, we consider the jet of higher Et- The jets are defined by the fc^-algorithm with the 
Snowmass merging rule. 

We shall also study the cross sections da /dx pn ', the variables x p ' 7 being either the x b s 
variables 

P't 



oos 2EP'"/ 

or the xll variables 



A P ±v h _l et P ±v jet 



n h (p ± v h + p ±n jet 
2Ep>~( 

E p ^ are the energies of the incoming proton respectively photon and the plus sign in 
e ±rj corresponds to x p , the minus sign to x 1 . The main difference between x b s and xll 
consists in the fact that the definition of xll does not require the measurement of the jet 
transverse energy. Further it has to be noted that - contrarily to the dijet cross section 
- the variables defined in (|6|) and ([?]) do not even for the Born contribution coincide with 
the variables Xp arton and x^ arton in the parton densities of the photon and the proton. The 
latter are (for the Born contribution) given by 

h ( ± v h + e ±v ]et ) 

■^parton 2E P '"< ' 

Therefore the partonic variables x^ rton are larger by a factor 1/z than the variables 
and larger by a factor (e ±T,h + e ±vjet )/(z e ±r]h + e ±7?3e *) than the variables x^ s . At NLO, 
where a third parton in the final state is involved, the situation is more complicated. 
The variables xj^J and x v £l defined in (g) and (|?j) differ by about x^ s /x v £l ~ (z + l)/2z 
at central rapidities. This difference will be discussed below in the context of Figs. || and 

Fig. |6| displays the direct and resolved contributions, calculated with the scales \x = 
M = Mp = Pt, as a function of 7] h , integrated over \rf et \ < 2. The two contributions are 
comparable, except in the forward region where the resolved one is much larger. However, 
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Figure 6: Rapidity distribution for the hadron-jet cross section da/drf 1 at the scales fi 
M = Mp =p T , integrated over 7GeV< p T <15GeV, El et > 5 GeV, \^ et \ < 2. 



we must keep in mind that the pattern of the separate contributions to the cross section 
depends on the choice of the scales /i, M and Mp ; only the sum of the resolved and direct 
contributions has a physical meaning and can be compared to data. 

In Fig. [7] we display the scale dependence of the cross section dajdrf 1 by varying 
all scales simultaneously, \x = M = Mp = Cpp, with the parameter C in the range 
1/2 < C < 2. We see variations which are similar to those observed in Fig. |2| in the 
inclusive case for pp ~ 7 GeV, where the cross section varies between +30% and -20% for 
1/2 < C < 2. However, varying all three scales in the same way is a very rough measure of 
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Figure 7: Scale dependence of da/drj h (|ry je *| < 2). 

the scale dependence. A better study has been given in Table § where we showed that the 
cross section exhibits a plateau in the scales /i and M for 0.3< <0.5 and 0.5< Cm <2. 
We performed a similar investigation for the hadron-jet cross section with the kinematic 
conditions of Fig. [7|. However, contrarily to the preceding study, we look for an optimum 
of the cross section in all three scales fi, M and Mp. That is, we look for optimal scales 



where dcr/d/j, = da/dM = da/dMp = \21\. We performed this complete study for 



da jdrf 1 with rj h in the range 0.5 < rf 1 < 1. The results are displayed in Table ||. 

This investigation is very CPU time consuming. Therefore, for a study of the optimum 
of do jdrf 1 in various bins in rj h , we used the following simplified approach. We fix the scale 
Mp to the optimal value of Table |3|, Mp = 0.5 pp, and look for optima of da/drf 1 in the 
(C^, CM)~plane for the rf 1 bins of Fig. |7j. The corresponding cross section is shown in Fig. |7] 
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( J 71 /f 






du j dff 1 [nb] 


1 


~ 0.21 


~ 1 


0.182 


0.5 


~ 0.3 


~ 1.5 


0.173 


0.3 


~ 0.3 


~ 1.5 


0.180 



Table 3: Scale optimisation for the hadron-jet cross section integrated over 0.5 < r/ 1 < 1, 
\i]i et \ < 2, 7GeV < p T < 15GeV and E J T et > 5GeV. 

and we see that it is close to the cross section obtained with p, = M = Mp = px/2. It is 
encouraging to note that the data of Fig. |l] are in good agreement with theory calculated 
with C = 0.5 in the range pr > 7GeV. 

Let us now discuss the cross sections da/dx b s and da/dxLL- Features of this cross 
sections are displayed in Figs. || to |l4[ Fig. ^ gives a clear illustration of the differences 
between the variables da/dx b s and da/dxn- In the partonic variable Xp arton (see eq. ®), 
the cross section would have a peak at Xp ar -f- on ~ 1, mainly due to the direct Born con- 
tribution which is proportional to 5(1 — Xp arton ). This peak is shifted to lower values of 
x\ L and xl bs due to the relation p\ = zp^, where a is the "parent" parton of the hadron 
h. This shift is larger for the variable xl L , as discussed at the beginning of this section. 
We observe a similar, but much less pronounced pattern for the proton variable x p ; the 
partonic distributions in the proton decrease rapidly with x v varton and this behaviour is 
reflected in the distributions of Fig. ||, da jdx\ L having a smaller width than da/dx p ohs . 
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Figure 8: Comparison of da/dx2 bs and da/dxJ^ L , integrated over 7GeV< <15GeV, 
E 3 T et > 5GeV, \rf et > h \ < 2. 
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Figure 9: Comparison of da/dx^ bs and da/dx p LL) integrated over 7GeV< px <15GeV, 
E T et > 5GeV, \rr> et ' h \ < 2. 

Fig. ^ displays a detailed study of the resolved and direct contributions as functions 
of x\ L and x2 bs - As already remarked, the shift of the peak to lower values of x 7 is 
more pronounced for xJ jL than for xL . The effect of cuts in 77" and t^ 6 * are also clearly 
visible, the kinematic region < rj ,rf et < 2 contributing to the low x 7 domain, and 
the region —2 < r] h ,rf et < to the large x 1 domain. The relative size of the resolved 
and direct contributions as given in Fig. is of course dependent on the scale choice 
(/x = M = Mp = pt for Fig. |To| ). As already said at the beginning of this section, only 
the total cross section is a physical observable. The variation of the size of the direct 



and resolved contributions with the scales varied by the factor C is shown in Fig. 11. For 
instance, we observe that for 0.5 < C < 1 the resolved component is almost stable whereas 
the direct component strongly varies. 
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Figure 10: Comparison of resolved and direct contributions to da jdx^ ohs and do jdx\ L in 
different rapidity ranges. 
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Figure 11: Comparison of resolved and direct parts for different scale choices, integrated 
over the rapidity range — 2 < r\ h ^ et < 2. 
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The hadron-jet cross section offers the possibility to measure the parton distributions 
in the proton and in the photon. The quark distributions are constrained by DIS exper- 
iments (very well in the proton case, and with rather large errors in the photon case). 
Therefore, concerning the quark distributions, the hadron-jet cross section can only put 
some additional constraints on F a / 7 (x, M). The situation with respect to the gluon dis- 




Figure 12: Contribution of the gluon in the photon to da/dx^ bs and dajdx\ L . The 
rapidities are integrated over the range — 2 < r] h ^ et < 2. 

tribution is different. The gluon distributions are not well determined in DIS experiments 
because the virtual photon couples at leading order only to quarks, such that the gluon 
distributions appear only at the level of higher order corrections. This is not the case 
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for the hadron-jet cross section where the gluon distributions appear already at the Born 
level. Therefore we expect an important sensitivity of the cross section to the gluon distri- 
butions, especially in the kinematic regions where x par ton is small. This point is illustrated 
in Figs. H to |J. 
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Figure 13: dcr/dxlr with the hadron and jet rapidity cuts 1 < r] fl ,rf et < 3. 



In Fig. 12, we observe that the gluon in the photon makes a large contribution at small 
x 7 , x 7 < 0.25. However, the cross section is small in this x 7 region because the cuts r] h < 2 
and r7 je * < 2 forbid to reach low values of x 7 . Using different cuts which reinforce the 
small x 7 region, namely 1 < r] h ,rf et < 3, we obtain the results shown in Fig. 13. We see 



that the cross section in the small x 7 region is much larger and made up almost entirely 
by the resolved contribution. The AFG02 parton distributions for the photon allow to 
modify the normalisation of the non-perturbative VDM component. In order to exhibit 
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the sensitivity to this component, we show a curve where the coefficient of this VDM input 
has been set to zero. We also show the magnitude of the gluon contribution to do jdx\ L 
in this kinematic range. We conclude from Fig. 13 that the rapidity cuts 1 < rj h ,rp et < 3 
select a kinematic region where the sensitivity of the hadron-jet cross section to the gluon 
distribution in the photon is very large. 

The contribution of the gluon in the proton turns out to be large even with rapidities 
integrated in the whole range —2 < r) h, ^ et < 2, as can be seen from Fig. [14]. We also 
observe that the effect of using different parton distributions for the proton, namely the 



MRST99 |22J, the new MRST01 []2J and the CTEQ6M |24| sets, is quite large and mainly 
due to the different shape of the gluon in the different sets. Although the variation of the 
cross section when varying all scales simultaneously between 0.3 < C < 2 is larger than 
the variation due to the different pdf sets, the scale variations rather produce an overall 
shift of the curve, but do not change the shape. On the other hand, the fact that the gluon 
distribution of MRST01 peaks at higher x values than the one of CTEQ6M (see Fig. 16 of 



[2"4|| ) is clearly reflected in Fig. 14. We also note that the region x\ L ~ 0.02, corresponding 
to 

x ^>arton ~ 0.05, is an interesting x-range because it lies in the window between constraints 
from HERA for lower x and the Tevatron jet data for higher x. Therefore the cross section 
do jdx\ L can serve to further constrain the parton distributions of the proton. 
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Figure 14: da/dx\ L calculated with different pdf sets and size of the gluon in the proton : 
rapidities integrated over the range — 2 < r\ h ^ et < 2. 
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5 Conclusion 



We have studied the photoproduction of inclusive large-pT charged hadrons and the pro- 
duction of a charged hadron plus a jet. For the inclusive case we compared the pr- and 
rapidity distributions to HI data |Kj and found reasonable agreement. However, for a 
value of p™ m as low as 3 GeV, the dependence of the NLO result on the renormalisation 
and factorisation scales is very large. Only for pji m >7 GeV a plateau where the cross 
section is approximately stable against scale variations could be found. We also studied 
the effect of using different fragmentation function parametrisations J?], H ^] and compared 
to a previous analysis of Kniehl, Kramer and Potter ||. For the parton distributions in 
the photon, we used the new AFG02 parametrisations [O]. 

For the hadron-jet cross section, we studied the rapidity distributions and the cross 
sections da / dx p , da / dx 1 . We analysed the difference between the partonic momentum 
fractions x^ rton and the observables x^ s defined via the observed transverse momenta 
and rapidities of the hadron and the jet. We further proposed a variable which does 
not require the measurement of the jet transverse energy. 

We also carried out an exhaustive study of the scale dependence. We found a stability 
region for the cross section integrated over 7 GeV < p^ < 15 GeV, E^ 1 > 5 GeV and 
performed a scale optimisation. 

Finally, we investigated the possibility to constrain the parton distributions (in par- 
ticular the gluon distributions) in the photon and in the proton via the hadron-jet cross 
section. We show how rapidity cuts can increase the sensitivity to the gluon distributions 
in the photon. We also found a rather large sensitivity to the parton distributions in the 



proton. We show a comparison of the MRST99 [22] distributions to the new MRST01 



and CTEQ6M [24] distributions. 
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